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1 Introduction 



The 77 (or 'di-photon') inclusive signature is possibly the most studied one experimentally, in 
the context of the Large Hadron Collider (LHC), since it allows for the detection of a relatively 
light ( < 130 GeV) neutral Higgs boson, both within the Standard Model (SM) and the Minimal 
Supersymmetric Standard Model (MSSM) [1, 2]. In fact, the existence of such a light particle state 
has been hinted already by LEP2 data, both as a possible resonance in the region MHi ggs ~ 115 
GeV and as the mass range that best accommodates the higher order Higgs boson contributions 
to precision Electroweak (EW) data (see, e.g., [3], for a review). 

Despite having a Branching Ratio (BR) that is only at the permille level, the di-photon 
signature is preferred to the one associated with the main Higgs decay channel in the above 
mass region, i.e., Higgs — > bb (into pair of 6-quark jets, with BR practically one), as the latter 
is swamped by the huge QCD background typical of a hadronic machine, whereas the former 
is much cleaner in such an environment. Besides, the higher precision that one can achieve in 
determining both directions and energies of the photons (as compared to those of jets), allows one 
to obtain a Higgs mass resolution of the order 2-3 GeV, which compares rather well to a typical 
15-20 GeV accuracy from jet events, also recalling that the Higgs width is at most a few tens of 
MeV in the above mass region (so that, the worse the mass resolution the larger the background, 
while the signal size remains relatively stable) [4]. 

The di-photon signature may also represent a distinctive signature of broken Supersymmetry 
(SUSY), namely, in Gauge Mediated Supersymmetry Breaking (GMSB) models [5]. In these 
scenarios, the Lightest SUSY Particle (LSP) is the so-called gravitino, G. Moreover, the next-to- 
LSPs (NLSPs) are usually the lightest neutralino, x{-> or the lighter stau, ti, depending upon the 
actual configuration of the GMSB model. In the first case then, the following decay chain could 
well occur: Higgs — ► XiX^j ~ * 77 + ^miss = 1,2), where the missing (transverse) energy is 
due to the two LSPs and possibly neutrinos (from \2 —> Xi + vv) escaping detection 4 . 

Current collider limits however forbid the lightest MSSM Higgs boson, h, to decay into two 
neutralinos, i.e., Mh < M~ +M~ (i,j = 1, 2), so that only the two heavier neutral Higgs bosons, 

Xi Xj 

H,A, can initiate the above SUSY decay chain. Recalling that h —> 77 direct decays can still 

occur in GMSB models, it is intriguing to consider the possibility that all neutral Higgs bosons 

of the MSSM can be detected in the same channel, that is, a pair of photons accompanied by 

4 Notice that, within the MSSM, one can obtain the di-photon signature (including -E m i ss ) of heavy Higgs bosons 
also from non-GMSB scenarios in which the LSP is the lightest neutralino, X?, such as in (minimal) Supergravity 
(mSUGRA) models. In fact, the following decay chain can occur here: H/A — > X2X2 followed by a double radiative 
decay X2 - * X?7 [6]- It has been shown in Ref. [7] that this channel can be large if tan/3 (see below for its definition) 
is small (~ 1), with or without gaugino mass unification at high scales. However, this MSSM configuration is ruled 
out by the LEP2 limit on the lighter chargino mass [8]. Only under the assumption of non-universal gaugino masses 
can the BR(x1> — > Xi7) be- large at high tan (3 in mSUGRA models [9] . 
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some amount of missing (transverse) energy, E m i SS . Notice that the latter should in average be 
larger for signal events as compared to the background, which is dominated by prompt di-photon 
production, where E m i ss is mainly due to jet energy losses down the beam pipe or to non-fully 
hermetic detectors. In contrast, one would naturally expect a large E m i ss value arising from 
the above H and A decays (if not from invisible decay products of heavy particles produced in 
association with any of the Higgs states, see below), so that the missing (transverse) energy may 
be used in the kinematical selection. In practice though, only two resonances could in the end be 
visible, as the H and A bosons are almost degenerate in mass. Besides, the latter would tend to 
be located at higher 77 invariant masses as compared to the one due to h — > 77 decays, where 
the production cross section is smaller, but so is the di-photon continuum. 

Finally, whereas the h — > 77 resonance can directly be reconstructed from the photon four- 
momenta, the same is not true in the H,A — > XiXj ~^ 77 + E m - lss = 1,2) channel. Here, 
however, after ensuring that the two photons are not back-to-back, one can attempt to resolve the 
.E m i ss along their directions and add it to the photonic transverse momenta, and . Scaling 
up the respective sums by the ratios V 11 IVr an d V 12 IVr gi ves i n principle the reconstructed 
momenta of the neutralino pair. In practice though, the presence of several unresolved sources of 
missing (transverse) energy may spoil the mass reconstruction. 

It is the purpose of this paper to investigate in detail such a phenomenology, in the context 
of GMSB scenarios. After a brief discussion of the parameters defining GMSB models and a 
description of the tools used in order to carry out our numerical studies, we will present the 
results and draw our conclusions. 

2 The spectrum in GMSB models 

In GMSB models, the symmetry of the Superpotential is broken at some relatively low scale, say, 
a few hundred TeV (the 'hidden sector'), and SUSY-breaking is communicated to the detectable 
particles (the 'visible sector') through so-called 'messenger' fields, effectively, gauge bosons. 

In fact, renormalisability of the theory, coupled with economy of field content, dictates that the 
messenger sector (MS) be comprised of chiral Superfields such that their SM gauge couplings are 
vectorial in nature. Most GMSB models actually consider these fields to be in (5 + 5) or (10 + 10) 
representations of SU(5). They are also chosen to transform as a multiplet of a Grand Unification 
Theory (GUT), so that the SUSY prediction of gauge coupling unification is preserved. These 
facts restrict the maximum number of messenger families Nm to be < 4 and < 1 for the (5 + 5) 
and (10 + 10) constructs, respectively. 

Limiting ourselves, for the time being, to a single pair of MS Supermultiplets ( X I' + 1 5'), consider 
a term in the Superpotential of the form A5^^, where S is a SM singlet. The scalar (S) and 
auxiliary (F5) components of S may acquire Vacuum Expectation Values (VEVs) through their 
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interactions with the hidden sector fields. SUSY-breaking is thus communicated to the MS, with 
the fermions and sfermions acquiring different masses. This, in turn, is communicated to the 
SM fields resulting in the gauginos and sfermions acquiring masses at the one-loop and two-loop 
levels, respectively. The expressions, in the general case of multiple messenger pairs and/or gauge 
singlets Si, is a somewhat complicated function [5] of M = (S) and A = (Fs)/(S). However, if 
there is just one such singlet, the expressions for soft SUSY-breaking gaugino and scalar masses 
at the messenger scale M simplify to: 

M,(M) = Wm ^ 9 7*V (!) 



4vr 3 \MJ 

n>)(M) = 2iV„A^(A)g KiC /(^) 2 . (2) 

In eq. (2), C( are the quadratic Casimirs for the sfermion in question. The factors ki equal 1, 1 
and 5/3 for SU(3), SU(2) and U(l), respectively, with the gauge couplings so normalised that all 
HiCe^s are equal at the messenger scale. The threshold functions are given by 

1 + x 

gi(x) = — 2~log(l + x ) + ( x ~ x )i ( 3 ) 



92( X ) = o 



log(l+x)+2Li 2 (- 



x \ 1 / 2x 
Tx~) ~ 2 l2 \l + x 



+ (x -> -x). (4) 



The Superparticle masses at the EW scale are obtained from those in eqs. (l)-(2) by evolving 
the appropriate Renormalisation Group Equations (RGEs). For the scalar masses, the D-terms 
need to be added too. The Higgs sector of the 'minimal' GMSB model contains the two usual 
Higgs doublets of the MSSM, (H u ,Hd). The ratio of the VEVs of the latter is parameterised as 
tan/3 = Moreover, in the Superpotential one has a Higgs bilinear term of the form BuH u Hd- 
In general, fi and B depend on the details of the SUSY-breaking mechanism. However, if we 
assume that the EW symmetry is broken radiatively, then the values of /x 2 and B are determined 
in terms of the other parameters of the model. Without loss of generality, one may express the 
entire particle spectrum of such a GMSB model in terms of five external inputs only: M, A, tan (3, 
sgn(/x) and Nm (hereafter, we assume Nm = !)• 



3 Parameter scans 

As a starting point of our investigation, we shall discuss the relevant particle spectrum by searching 
for regions of the GMSB parameter space where the heavy Higgs bosons, H and A, and the lightest 
neutralino, Xii have masses and compositions (in terms of the gaugino-Higgsino fields) such that 
the decays H,A^ XiXi are kinematically allowed and reach BRs that are at least comparable 
to that of the SM Higgs decay into two photons, for which one has BR(i?sM — ► 77) — 10~ 3 . 
In fact, such a channel is of extreme importance in the MSSM too, as already noted. Here, the 
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parameters tan (3 and Ma define entirely the Higgs sector of the MSSM at tree-level, tan (3 having 
being already defined and with Ma being the mass of the pseudoscalar Higgs state A (the other 
two states, h and H, are scalars) 5 . It turns out that, for any tan/3 value, if Ma ~ 150-200 GeV, 
the h — > 77 decay mode is a discovery channel of the lightest Higgs boson of the MSSM. 

In order to sample the strength of the di-photon signal of our interest, we present a set of 
numerical results that include the Higgs and lightest neutralino masses, as well as the BRs of the 
channels H, A — ► XiXi- In computing these rates we have used the subroutines of IS A JET 7.58 
[10] that implement the GMSB model, with several choices of parameter inputs. Namely, in Tabs. 
1,2 and 3 we have fixed M = 150 TeV and taken A = 75, 100 and 125 TeV, respectively. Here, 
the sign of fi is always negative whereas tan (3 varies from 5 up to a maximum value where the 
decays H,A^ XiXi are no longer kinematically allowed. 

Prom Tabs. 1-3 one can appreciate the following trends. 

• As intimated in the Introduction, the Higgs masses M# and Ma show a degeneracy within 
a couple of GeVs. Thus, we can add their corresponding two-photon signals, as it would 
not be possible to distinguish among them solely on the basis of the reconstructed mass. 

• The lightest neutralino Xi decays into a photon plus a gravitino, xi ~^ 7 + G, while 
BR(H, A — > XiXi) can exceed the 10~ 3 level, this yielding altogether a decay rate for 
the channels H,A^ XiXi ~~ * 77 + -E'miss above our reference di-photon SM Higgs decay 
rate. 

• For values of tan (3 larger than 40 or so, the decays H,A^ XiXi are no longer kinematically 
allowed. 

Table 1: Higgs, lightest neutralino masses and BR(H, A — > XiXi) f° r a sample set of GMSB inputs 
with (M, A) = (150,75) TeV. 



tan [3 


M H [GeV] 


M A [GeV] 


M~o [GeV] 


BR(H) 


BR(A) 


5 


430 


428 


105 


1.4 x 10~ 2 


1.1 x 10~ 2 


10 


404 


403 


104 


1.1 x 10~ 2 


1.2 x 10- 2 


15 


392 


391 


103 


6.0 x 10~ 3 


7.3 x 10- 3 


20 


377 


377 


103 


3.5 x 10~ 3 


4.6 x 10~ 3 


25 


358 


358 


103 


2.1 x HT 3 


3.0 x 10~ 3 


30 


332 


332 


103 


1.3 x HT 3 


2.0 x 10~ 3 


35 


296 


296 


103 


7.5 x 10" 4 


1.4 x 10~ 3 


40 


245 


245 


102 


2.5 x 10~ 4 


7.9 x 10~ 4 


45 


162 


161 


102 









5 In fact, in the GMSB model, all Higgs masses are derived quantities, as specified in the previous Section. 
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Table 2: Higgs, lightest neutralino masses and BR(H, A — > X1X1) f° r a sample set of GMSB inputs 
with (M, A) = (150, 100) TeV. 



tan/3 


M H [GeV] 


M A [GeV] 


M~o [GeV] 


BR(H) 


BR(A) 


5 


552 


551 


146 


7.6 x 10" 3 


7.5 x 10" 3 


10 


521 


520 


145 


7.1 x 10" 3 


8.6 x 10" 3 


15 


505 


505 


145 


3.7 x 10" 3 


4.9 x 10" 3 


20 


488 


487 


144 


2.1 x 10" 3 


3.0 x 10" 3 


25 


464 


464 


144 


1.2 x HT 3 


1.9 x 10~ 3 


30 


433 


432 


144 


7.6 x 10" 4 


1.3 x HT 3 


35 


390 


389 


144 


4.1 x 10~ 4 


8.6 x 10~ 4 


40 


330 


330 


144 


1.2 x HT 4 


4.7 x 10~ 4 


45 


237 


237 


144 









Table 3: Higgs, lightest neutralino masses and BR(H, A — > XiXi) f° r a sample set of GMSB inputs 
with (M, A) = (150, 125) TeV. 



tan p 


M H [GeV] 


M A [GeV] 


M~o [GeV] 


BR(H) 


BR(A) 


5 


666 


664 


197 


4.8 x 10" 3 


5.4 x 10" 3 


10 


628 


628 


195 


4.8 x 10" 3 


6.5 x 10" 3 


15 


610 


610 


195 


2.5 x 10" 3 


3.7 x 10" 3 


20 


590 


589 


195 


1.4 x 10" 3 


2.2 x 10" 3 


25 


562 


562 


195 


7.9 x 10" 4 


1.4 x 10" 3 


30 


526 


526 


195 


4.4 x 10~ 4 


9.3 x 10~ 4 


35 


477 


476 


195 


2.1 x 10~ 4 


5.9 x 10~ 4 


40 


408 


408 


194 


2.3 x 10~ 5 


2.3 x 10~ 4 


45 


306 


305 


194 









For some regions of the GMSB parameter space, the decays H,A — > XiX^X^X^ followed 
by X2 ~~ *■ Xi + X can be significant enough to contribute to the two-photon signals. Although 
this effect was not included in the previous Tables, as it was small, it will be considered in the 
remainder of this Section and in the Monte Carlo (MC) simulations of the next one as well. To 
discuss their effects, we have produced three additional sample points, which we elaborate upon 
below and in Tabs. 4-6, before moving on to the event generator analysis. 

Tab. 4 corresponds to the Snowmass slope M = 2A [11], as we have fixed tan/3 = 15 and 
taken sgn(^i) positive, further varying A from 75 to 150 TeV (for lower values of A we get a 
chargino lighter than 150 GeV, which is not allowed by current collider bounds [12]). Both the 
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BR columns (BR(i7) and BR(A)) contain three rows, each corresponding from top to bottom to 
the decay rates of H and A into x?X?> X1X2 an d X2X2 pairs, respectively. 

For these sets of GMSB parameter space points the BR of H and A into X1X2 always dominates 
over the other two channels. This can be understood in terms of the enhancement of the H,A — 
Xi ~ X2 coupling, which over-compensates the phase-space suppression in the H,A^ X1X2 decay 
modes. Furthermore, one can see from Tab. 4 that, as A increases, the Higgs masses Mh and 
Ma become rather heavy and they still show degeneracy. The discussed BRs stay above 10~ 3 , 
but because the production rate decreases substantially for Higgs masses above about 600 GeV, 
only the lower values of A (~ 75 TeV) will produce sizable signals. It is interesting to note that 
in this case the BR of the lightest Higgs boson into two photons remains close to the SM value, 
with M h ~ 115 GeV. 



Table 4: Higgs, lightest and second lightest neutralino masses and BR(H,A —> x?Xi> X1X21 X2X2) 
for a sample set of GMSB inputs with M = 2A, tan (3 = 15 and sgn(jit) positive. 



A [TeV] 


Mu [GeVl 


Ma \GeV] 

■ LV± A. v J 


M~o \GeV\ 


M~n [GeVl 


BR(H) 


BR(A) 


75 


395 


394 


101 


184 


7.5 x 1(T 3 
2.08 x 10~ 2 
4.01 x 10 -3 


1.0 x 10~ 2 
4.4 x 10~ 2 
3.3 x 10~ 2 


80 


419 


418 


108 


198 


6.5 x 1(T 3 
1.8 x 1(T 2 
3.05 x 10- 3 


9.3 x 10~ 3 

4.1 x 10~ 2 

3.02 x 10~ 2 


90 


467 


466 


123 


225 


5.2 x 1(T 3 
1.56 x lO" 2 
1.56 x 10 -3 


7.6 x 10~ 3 

3.7 x 10~ 2 
2.4 x 10~ 2 


100 


515 


514 


137 


252 


4.2 x 10~ 3 
1.32 x 10~ 2 
6.47 x 10" 4 


6.2 x 10~ 3 
3.25 x 10~ 2 
1.74 x 10~ 2 


125 


631 


631 


173 


320 


2.8 x 10~ 3 
9.2 x 1(T 3 


4.3 x 10~ 3 
2.55 x 10~ 2 


150 


745 


744 


209 


387 


2.0 x 1(T 3 
6.9 x 10~ 3 


3.2 x 10~ 3 
2.1 x 10~ 2 
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Table 5: Higgs, lightest and second lightest neutralino masses and BR(H, A —> XiXii XiX2> X2X2) 
for a sample set of GMSB inputs with M = 2A, A = 75 TeV, tan (3 = 10 and sgn(/x) positive. 



A [TeV] 


M# [GeV] 


M A [GeV] 


M~o [GeV] 


M~o [GeV] 


BR(iT) 


BR{A) 


75 


405 


404 


101 


183 


1.54 x 10~ 2 
4.36 x lO- 2 
1.2 x 10^2 


2.0 x 10^2 
8.5 x 10-2 
7.0 x ID" 2 



Table 6: Higgs, lightest and second lightest neutralino masses and BR(H,A — > XiXi; XiX2> X2X2) 
for a sample set of GMSB inputs with M = 2A, A = 75 TeV, tan j3 = 35 and sgn(^) negative. 



A [TeV] 


M H [GeV] 


M A [GeV] 


M~o [GeV] 


M~o [GeV] 


BR(iJ) 


BR{A) 


75 


296 


296 


102 


188 


7.53 x 10~ 4 
9.73 x 10- 5 


1.4 x 10~ 3 
2.0 x 10~ 3 



4 Monte Carlo simulation of signal and background 



In this section, we evaluate the inclusive production cross sections at the LHC for the two heavy 
neutral Higgs bosons of the MSSM, H and A, each followed by all possible decays yielding two 
photons and missing (transverse) energy. 

Hereafter, we make the assumption that coloured SUSY particles (chiefly, squarks) are heavy 
enough 6 so that they do not enter the loops in the l gg — > Higgs' production mode nor they can 
produce Higgs bosons in cascade decays or enter the Higgs decay chains (the same for gluinos). 
Alongside the mentioned gluon-fusion channel, we consider the following Higgs production modes: 
'qQi99 QQ Higgs' (associated heavy-quark production), l qq — > qqV*V* — > qq Higgs' (vector- 
fusion) and 'qq( — ► V Higgs' (Higgs-strahlung) , with q representing all possible combinations 
of light (anti)quarks and where Q = b, t, V = W^,Z and Higgs = H,A (except in the last 
two modes, where the pseudoscalar Higgs boson cannot be produced). These are the leading 
production modes of neutral Higgs states of the MSSM at the LHC (under the above assumption 
of heavy squarks and gluinos) 7 and have been computed here at next-to-leading order accuracy, 
by adopting the programs described at [14], with default settings. 

As for the calculation of the SUSY decays rates (again, produced with ISAJET v7.58), we 
have adopted the GMSB configurations given in Tab. 7 for M and A, further varying tan (3 from 
5 to 40, always with sgn(/x) > 0. (Note that all the above choices of SUSY parameters are allowed 
by the lighter chargino mass limit, as derived within the GMSB model.) 

6 As it is typical in most GMSB scenarios. 
7 See Ref. [13] for a review of their properties. 
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Table 7: Sets of GMSB parameter points M and A considered in Fig. 1. 



Point 


M [lev J 


A [levj 


A 


150 


^C 

ID 


i ) 

D 


150 


100 




i 

10U 


IOC 

I/O 


D 


200 


75 


E 


200 


100 


F 


200 


125 


G 


200 


150 


H 


200 


175 



We define the 'effective' production rate of the 77 + £miss signature, a c g, as: 

o-cff = a(gg ^ H)x BR eS (H) 

+ a(gg -> A) x BR cfr (A) 

+ a{qq,gg^QQH) xBR cS (H) 

+ ^(99,55 -»• QQA) x BR eff (A) 

+ a(qq^HV) X BR eff (#) 
+ 99-H") x BR cS (H) 

where BR e g(H) and BR e g(^4) are defined as follows: 
BR eS (H/A) = 



(5) 



BR(H/A -> X1X1) + BR^M ^ Xixl) x Br(x° - x? + £ 

+Br(#M - X2X2) x BR.(x° - Xi + £ w) 2 ] x BR(x? - 7G) 2 .(6) 

In Fig. 1 we display the variation of a c s with tan (5 for the sets of GMSB parameter space points 
given in Tab. 7. From the pattern of each curve it is clear that the lower A the higher a e g. 
Furthermore, at low tan (3, a c s is larger than at high tan (3. This can be understood from the 
fact that, as tan (3 grows, the lighter stau (77) becomes the NLSP, with the decay Xi ~^ + r=F 
becoming the dominant one: this explains the sharp fall in cr c fj at large tan (3. 

The event simulation has been carried out by exploiting the SUSY implementation [15] of the 
HERWIG event generator [16], supplemented by the ISASUSY [10] subroutines (v7.58) interfaced 
to the ISAWIG code [17] for SUSY spectra generation. We list here the series of process numbers 
used for the MC event generation: i.e., 

IPRDC = 3320 3360 3375 3630 3815 3826 
3310 3325 3365 3610 3710 3816 3835 
3315 3335 3370 3620 3720 3825 3836 
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for the signal and IPROC = 2200 for the background (corresponding to direct di-photon production 

qq,gg -> tt) 8 - 

As illustrative values for the MC simulation, we have used the three points given in Tab. 8. 
Notice that they have already been discussed, as they are those appearing in: the fourth row of 
Tab. 4 (Point 1, which is SPS8 of Ref. [11]), Tab. 5 (Point 2) and Tab. 6 (Point 3). 

Table 8: Sets of GMSB parameter points used in the HERWIG simulation and to produce 
Figs. 2-4. 



Point 


M [TeV] 


A [TeV] 


N M 


tan p 


sgn(/x) 


1 


200 


100 


1 


15 


+ 


2 


150 


75 


1 


10 


+ 


3 


150 


75 


1 


35 





The signature we are trying to extract is simply defined as follows, along the lines of the 
ATLAS/CMS triggers [1, 2]. 

• Two photons are required: one with px > 40 GeV and 20 GeV for the other, both within 
2.5 in pseudorapidity. 

• The two photons are required to have a relative angle less than 175 degrees, in order to 
enable the mass reconstruction of the two heaviest Higgs bosons decaying into two photons 
and two LSPs, the latter yielding the missing (transverse) energy. 

• No cuts in missing (transverse) energy are enforced at this preliminary stage, nor any re- 
striction on the underlying hadronic activity is imposed. 

We then look at: 

• -Emiss: the missing (transverse) energy. 

This quantity is plotted in Fig. 2(a,b,c), normalised to one. (Hereafter, the labelling a(b)[c] in the 

figures refers to GMSB set 1(2) [3] as given above in Tab. 8.) As intimated in the Introduction, one 

may appreciate here the fact that the missing (transverse) energy distribution is much softer for 

the background, as compared to the signal. A suitable cut on this quantity, which will enhance the 

8 The attentive reader will notice that we have generated in the MC simulations more processes than those 
used in Fig. 1. This has been done for completeness mainly, as the four channels described in Sect. 3 are indeed 
those making up most of the visible cross section. Furthermore, all possible decay channels of Higgs bosons and 
neutralinos are included in the MC simulation, through the ISAWIG input files. Also notice that the (inclusive) rates 
in Fig. 1 use next-to-leading order normalisation, whereas those in Figs. 2-4 adopt the lowest order one, as default 
in a MC event generator. 
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signal-to-background ratio, could be, e.g., E m i ss > 20 GeV. This may penalise signal contributions 
due to h decays, yet it will have a beneficial impact in extracting those arising from H and A. 

Upon enforcing this constraint, in addition to the cuts in 1.-2., we look at two kinematic 
observables: 

• m 77 : the invariant mass of the photon pair obtained by using their visible momenta; 

• M 77 : the invariant mass of the photon pair obtained after resolving the E'miss along the 
visible photon directions. 

These are plotted in Figs. 3(a,b,c) and 4(a,b,c), respectively, normalised to the integrated cross 
section in picobarns, as given by HERWIG. From these plots it is clear the potential of the LHC in 
detecting di-photon signals of neutral Higgs bosons, as induced by our sample of GMSB model 
configurations. However, while the direct h — > 77 resonance is clearly visible and narrowly centred 
around Mh in both distributions m 77 and M 77 (and so it was even prior to the enforcement of 
the cut in E m i ss ), no peak associated to H/A — > 77 + £" m i ss decays appears, yet the corresponding 
events spread over the entire mass range well above the di-photon background. In fact, our 
mass reconstruction procedure fails because of the many unresolvable sources of missing energy 
appearing at hadron level, once all decay modes of each unstable particle are allowed, as dictated 
by the MSSM configuration induced by the underlying GMSB scenario. Here, the exploitation 
of a more exclusive final state would be helpful to the cause of extracting the H/A resonance. 
However, we refrain here from pursuing this matter further and simply remark that, even in 
absence of heavy Higgs mass reconstruction, a clear excess in the di-photon channel above the SM 
expectations should be established at the LHC after minimal luminosity 9 , with some dependence 
in both the m 77 and M 77 spectra upon the relevant masses M h /a, M~ and Mq, whose actual 
values may be investigated via a comparison of the data to the MC generated distributions. 

5 Conclusions 

In summary, we have proved that, for some rather natural configurations of the GMSB parameter 

space consistent with current collider limits, one may be able to extract di-photon signals of all 

neutral Higgs bosons of the MSSM at the LHC. While, after customary ATLAS/CMS cuts on the 

two photons, the mass of the lightest state would always be visible in the form of a 77 resonance, 

the presence of the two heaviest states (which are degenerate in mass) can be established, in 

the form of a clear excess in the total number of 77 + -Emiss events over the corresponding SM 

predictions, after an additional threshold in missing (transverse) energy is enforced. Thus, after 

9 Recall that higher order QCD corrections to the background are well under control in comparison to the excesses 
seen in the last two figures, as they are of order 10-20% [18]. Moreover, the contribution to the background due to 
<19 —> 97 with the final state light (anti)quark mistagged as a photon, not included here, is also small in comparison. 
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such a selection, in order to test very specific GMSB model predictions, it suffices to exploit the 
total event rate. Finally, the study of the kinematics of the entire event sample may allow for the 
determination of crucial sparticle masses, such as those of the the LSP and NLSP, hence enabling 
one to strongly constrain the underlying SUSY-breaking mechanism. 

Our conclusions are based on a sophisticated MC event simulation but a more rudimentary 
emulation of detector response (based on Gaussian smearing of the visible tracks). However, our 
preliminary results are rather encouraging and we do believe that they call for attention on the 
ATLAS/CMS side, as the 77 channel is possibly the most studied one in the context of Higgs 
boson searches at the upcoming CERN hadron collider. 
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Figure 1: Variation of the effective cross section of eq. (5) with tan [3 for the representative points 
in the GMSB parameter space given in Tab. 7. 
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Figure 2: Differential distributions in -E m i ss (as denned in the text) normalised to one, after the 
cuts 1.-2. described in Sect. 4, as obtained from GMSB set 1 (top), 2 (middle) and 3 (bottom) 
defined in Tab. 8. The dashed(dotted) line represents the signal(background) rates. Bins are 2 
GeV wide. 
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Figure 3: Differential distributions in m 77 (as defined in the text) normalised to pb, after the 
cuts 1.-2. described in Sect. 4 and -E m i ss > 20 GeV, as obtained from GMSB set 1 (top), 
2 (middle) and 3 (bottom) defined in Tab. 8. The dotted(solid) line represents the back- 
ground (signal+background) rates. Bins are 3 GeV wide. 
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Figure 4: Differential distributions in M 77 (as denned in the text) normalised to pb, after the 
cuts 1.-2. described in Sect. 4 and E m i ss > 20 GeV, as obtained from GMSB set 1 (top), 
2 (middle) and 3 (bottom) defined in Tab. 8. The dotted(solid) line represents the back- 
ground (signal+background) rates. Bins are 3 GeV wide. 



